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We propose a generalization of the discrete Klein-Gordon models free of the Peierls-Nabarro 
barrier derived in Nonlinearity 12, 1373 (1999) and Phys. Rev. E 72, 035602(R) (2005), such that 
they support not only kinks but a one-parameter set of exact static solutions. These solutions can 
be obtained iteratively from a two-point nonlinear map whose role is played by the discretized first 
integral of the static Klein-Gordon field, as suggested in J. Phys. A 38, 7617 (2005). We then discuss 
some discrete </> 4 models free of the Peierls-Nabarro barrier and identify for them the full space of 
available static solutions, including those derived recently in Phys. Rev. E 72 036605 (2005) but 
not limited to them. These findings are also relevant to standing wave solutions of discrete nonlinear 
Schrodinger models. We also study stability of the obtained solutions. As an interesting aside, we 
derive the list of solutions to the continuum </> 4 equation that fill the entire two-dimensional space 
of parameters obtained as the continuum limit of the corresponding space of the discrete models. 

PACS numbers: 05.45.-a, 05.45.Yv, 63.20.-e 
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INTRODUCTION AND SETUP 



Discrete nonlinear models play a very important role in 
many physical applications [3, 0| ■ An important class of 
these models consists of a few completely integrable lat- 
tice equations -such as the Toda lattice Q, the Ablowitz- 
Ladik lattice jj], and the integrable sine-Gordon lattice 
H . The fact that these lattices possess exact soliton so- 
lutions demonstrates that, in principle, discreteness of 
the host media does not preclude the propagation of lo- 
calized coherent structures. Moreover, the mobility of 
soliton-like excitations in discrete media is a key issue in 
many physical contexts; for instance, the mobility of dis- 
locations, a kind of topological solitons, is of paramount 
importance in the physics of plastic deformation of met- 
als and other crystalline bodies 

A prototypical class of discrete models, relevant to a 
variety of applications, consists of the Klein-Gordon dy- 
namical lattices Q. One of the main representatives of 
this family of models is the so-called cf> 4 model |7j , which 
features a cubic nonlinearity. This simple power law non- 
linearity renders this model a ripe testbed for studying 
the existence and stability of nonlinear solutions, and 
comparing their properties in continua and lattices. 

In the (1 + l)-dimensional continuum framework (and 
in the absence of spatially dependent external poten- 
tials), a solution can be shifted arbitrarily along x by 
any Xo (x is the spatial coordinate and Xq = const.), due 
to the existing translational invariance. On the other 
hand, in the discrete system, translational invariance is 
generically lost and equilibrium static solutions exist only 
for a discrete rather than for a continuum set of xo Q- 
Some of these equilibrium solutions correspond to energy 
maxima and arc unstable, while others, corresponding to 
energy minima, are stable. The difference between such 
maxima and minima of the energy is typically referred to 



as the Peierls-Nabarro barrier (PNb) . 

It is of particular interest to develop discretizations 
that do not feature such barriers. In such cases, one 
might expect that the ensuing models would be more 
faithful representations of their continuum counterparts, 
regarding both symmetry properties and travelling solu- 
tion features (even though there are some notes of caution 
that should be made; see e.g. the discussion of [8j). 

In that vein, recently, a number of non-integrable dis- 
crete Klein-Gordon equations free of the Peierls-Nabarro 
barrier (PNb-free) has been systematically constructed. 
The first set of such models which were, by construc- 
tion, Hamiltonian ones, was obtained by Spei ght and 
co-workers |9( using the Bogomol'nyi argument |10|. A 
second successful attempt led to the construction of 
momentum-conserving discretizations [Til ]. It was then 
demonstrated, surprisingly, that the PNb-free models of 
that kind conserving both energy and classical momen- 
tum do not exist [IJ. New PNb-free cj> 4 lattices were 
derived by Barashenkov et al. [fa]. However, we note 
that for the lattices derived in [ij ^| , only the kink- 
type solution has been considered, while it is well-known 
that the continuum Klein-Gordon equation can support 
a number of other solutions [T^| . 

On the other hand, a general approach to the construc- 
tion of the PNb-free lattices was recently reported in . 
This approach is based on the use of the discretized first 
integral (DFI) of the corresponding static field equation, 
and the integration constant that enters DFI generates 
the one-parameter solution space. The DFI is in fact 
a two-point nonlinear algebraic equation from which the 
exact static solutions of the three-point PNb-free discrete 
models can be found. 

In this work we systematically use the DFI approach 
and derive the Hamiltonian and non-Hamiltonian PNb- 
free Klein-Gordon lattices supporting the one-parameter 
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space of solutions, generalizing the lattices offered in 
We then focus on the cj> 4 field and discuss in detail 
the ener gy-c onserving PNb-free model proposed very re- 
cently in |16| , as well as the momentum-conserving model 
of [llj (the two are identical in the static case) . For these 
models we describe the full solution space of the underly- 
ing static problem that also includes the solutions derived 
earlier in |16(. 

Our setting is the following: We consider the Hamil- 
tonian of the Klein-Gordon field, H = Ex + Ep, with 
the kinetic and potential energy functionals respectively 
defined as, 

1 f°° 

Ek = ^J tfdx, (1) 
E P = - [<f>l + Wt4>)]dx, (2) 

* J — oo 

where (f>(x, t) is the unknown field and V(<fi) is a given po- 
tential function. The corresponding equation of motion 
is 

<j>u=^-V l {(t>)=D{(t>{x ] t)), (3) 

where V'{<fi) = dV/dcfi. 

Equation will be discretized on the lattice x — nh, 
where n = 0, ±1, ±2..., and h is the lattice spacing. 

We would like to construct a nearest-neighbor discrete 
analog to Eq. © of the form 

<f) n = D(h,(f>n-l,4>n,<f>n+l)i ( 4 ) 

such that in the continuum limit (h — » 0) we have 

D(h, (f> n ^\,(j) n ,(j) n+1 ) -> D{4>), (5) 

and that the solution to the three-point static problem 
corresponding to Eq. 

D(h, n _i,0 n , <f> n +i) = 0, (6) 

can be found from a reduced two-point problem 

L/(M«-i,<W =0. (7) 

If this reduction is achieved, then the exact static solu- 
tions can be constructed upon solving the algebraic Eq. 
|J7)| iteratively, starting from arbitrary admissible value 
of 4> n —i or <j> n . Arbitrariness in the choice of the initial 
condition implies that the static solution can be placed 
anywhere with respect to the lattice and, for that reason, 
such lattices are called translationally invariant. 

Discretization Eq. |@J may result in a non-Hamiltonian 
model due to the non-potential nature of the background 
forces. In the latter setting, the absence of an energy 
functional renders ambiguous the definition of the PN 
barrier, hence we clarify this point in what follows. Sup- 
pose that we have two equilibrium solutions, <fin^ and (jffl 



(often, in Hamiltonian models with PNb, the second one 
corresponds to a linearly unstable energy maximum while 
the first one to a linearly stable energy minimum). The 
work done by the inter-particle and background forces to 
move the nth particle from the position <f>n to the posi- 
tion (f>n ) is W n = J D(h,(t> n -i,4> n ,4> n+ i)<l(f> n and the 

total work performed to "transform" the solution <fi n to 
(fin^ is W — X^oo Wn- The PN barrier is defined to be 
equal to W . For Hamiltonian models this definition is 
equivalent to the classical one because W is equal to the 
energy difference between the second (possibly higher en- 
ergy unstable) configuration and the first (possibly lower 
energy stable) configuration. For non-Hamiltonian mod- 
els, W n will depend on the path connecting the initial 
and final configurations of particles. To calculate W one, 
therefore, has to specify such a path. While the height 
of the PN barrier in a non-Hamiltonian lattice depends 
on the path, the PNb-free, non-Hamiltonian lattice can 
be unambiguously defined as the translationally invariant 
lattice where the quasi-static transformation between the 
two configurations of interest can be done continuously 
through the set of equilibrium configurations. Along the 
path through the equilibrium configurations, forces act- 
ing on particles are zero and thus, W n — for any n 
which, in turn, results in W = 0. This notion of the PNb- 
free lattice is also applicable to the Hamiltonian lattices 
and will serve as our definition of the PNb-free models; 
notice, however, that we do not overlook the mathemat- 
ical subtleties involved in this definition, including the 
question of whether W is zero along all realizable paths 
connecting </4 and (fii? in the non-Hamiltonian mod- 
els, among equilibrium configurations, as is the case for 
their Hamiltonian siblings. These and related questions, 
including an appropriate modified definition of the rele- 
vant quantities for models such as those of the discrete 
nonlinear Schrodinger type, will be left for a future pub- 
lication. 

Our scope in what follows is to generalize the approach 
developed in ^jjj to show how to construct all possible ex- 
act static solutions of the models of [l5| and |lq (see also 
HUE!)- including those derived in |l6j. This will also 
lead us to introduce a number of solutions (both localized 
and extended ones) that, to the best of our knowledge, 
have not been discussed/ analyzed previously, such as the 
ones that will be termed "inverted" (see below). We will 
also discuss the stability of certain solutions among the 
obtained ones, for each of the models of interest (since 
their stability properties are different). 

We also note in passing that while our presentation 
will be geared towards the (fi 4 models, our results regard- 
ing the existence of solutions can equally well be applied 
to discrete equations of the nonlinear Schrodinger (NLS) 
type [l^, such as e.g. the Ablowitz-Ladik model. In 
particular, let us consider equations of the form: 

ilfin = To^n-l - 2V>n + lpn+l) + f{lpn, ~4>n), (8) 
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(where the overbar denotes complex conjugate) with 
f(ili n exp(i0),ip n exp(-i0)) = f(tp n , ip n ) exp(i0) and 
lim/i-,0 ffym VvO — ~ A|V'| 2 '0- Then, looking for stand- 
ing wave solutions of the form ip n — exp{i\t)<j> n , with 
4> n real, one would obtain the Klein-Gordon static prob- 
lem for the standing wave spatial profile <j) n . Hence all 
the discussion given below for the existence of such solu- 
tions can be appropriately translated in the existence of 
standing waves of the corresponding discrete NLS mod- 
els. The reader should be cautioned however that the 
stability properties in the latter context may differ [a rel- 
evant example will be discussed later in the text]. 

Our presentation will be structured as follows. In sec- 
tion[n] we develop the DFI approach and derive PNb-free 
Klein-Gordon models supporting one-parameter space of 
static solutions. In section ITTT1 we present some PNb-free 
(/> 4 lattices and describe some of their basic properties. In 
section IIVI we discuss the details of the construction of 
the general exact static solutions of the models both in 
the localized (hyperbolic function) and in the extended 
(general elliptic function) form, and for both signs of the 
nonlinearity prefactor A. A = 1 corresponds to the so- 
called defocusing case, while A — — 1 corresponds to the 
focusing case in the standard terminology of NLS equa- 
tions. Examples of the solutions are then given in section 
IVl In section IVII we analyze the stability of the ob- 
tained solutions. In section IYlII slow kink dynamics in 
the PNb-free models is compared numerically with that 
in the classical <fi 4 model. In section IV 1111 we give a com- 
plete list of bounded and unbounded exact solutions to 
the continuum cf> 4 field. In section llXl we summarize our 
findings and present our conclusions. 



and 



)(i) = ±V^+c-vWi = o, 



(13) 



of which only the first one will be discussed below. 

We then construct the DFIs corresponding to Eq. © 
and Eq. I|12fl . which are respectively given by 



U(h, (j) n --t,(f>n) 



-l) 2 



h 2 



-2^ n _ 1 ,0„)+C = O ) 



(14) 



v(h, 4> n -l,4>n) = ±- 



-V2V(0„_i,0„)-C = O, 



(15) 



where we demand that V(cj> n —i, <j) n ) — ► V(4>) in the con- 
tinuum limit (h — > 0). 



B. Momentum-conserving PNb-free lattice 

First we construct a PNb-free Klein-Gordon lattice us- 
ing the unchanged form of the first integral, Eq. Q , and 
corresponding DFI, Eq. I|14|) . Calculating dU/dx and 
multiplying the result by (dx/d<p)/2, we find 



1 dU 
YdS 



D(x). 



(16) 



II. PNB-FREE KLEIN-GORDON LATTICES 



Discretizing the left-hand side of Eq. (|16[) we obtain the 
lattice Klein-Gordon equation 



A. Discretized first integral 

Following the lines of the DFI approach of Ref. [l5j . 
we start from the first integral of the static Eq. ©, 



U{x) = ct>l-2V{(j>)+C = Q, 



(9) 



where C is the integration constant. The first integral 
can also be taken in modified forms |15|. e.g., as 

v{x)=p[g{(g)-g{2V-C)\ =0, (10) 



or as 



w(x)=p[g(<l>l + C)-g(2V)}=0, (11) 

where p and g are continuous functions and p(0) = 0. 
Note that v(x) and w(x) are equivalent only if g(£) = £. 
We will consider the case of p(() — C, g{£.) — i.e., the 
unchanged form of the first integral, Eq. ©, together 
with the case of p(Q — C, g{0 — VC, for which we have 
the two possibilities 



)(x) = ±<f> x - ^2V((f>)-C = , 



(12) 



U(h,(f> n ,(p n+ i) - U(h,4> n -i,4> n ) 



»n-l 



(17) 



whose static solutions can be found from the two-point 
problem, Eq. ((TH) . and thus, the lattice is PNb-free. 

The lattice of this type was first derived in ^l| where it 
was also demonstrated that it conserves the momentum 



(18) 



The integration constant C that appears in Eq. 114|) 
cancels out in Eq. (|1 7|) . This means that all kinds of 
static solutions derived from Eq. I|14l) for different C 
(kink solution corresponds to C = 0) will be the static 
solutions to one and the same Klein-Gordon lattice, Eq. 
(|T7|l . since it does not depend on C. 

The right-hand side of Eq. H17[) becomes non-singular 
when the discretization of the potential V(cf>) is a 
polynomial function having symmetry V((f) n -i, <p n ) = 
V(4> n ,4> n -i). The most general expression of this type 
was given in ^1 ■ 
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C. Energy-conserving PNb-free lattice 



D. Possible generalizations 



To construct Hamiltonian PNb-free lattices we dis- 
cretize not the equation of motion, Eq. but the 

Hamiltonian, Eq. (JTJ and Eq. J2J). We now use the 
modified first integral in the form of Eq. with the 
upper sign and rewrite the potential energy functional, 
Eq. ©, as follows, 



Ep — — 



v{x)} 2 +2^v^(Fc) dx, (19) 



where we omitted the constant term. Discrctizing the 
kinetic energy, Eq. (I), and the potential energy, Eq. 
()19p. we obtain the discrete Hamiltonian 



+ 2 



) - C} ■ (20) 



If the background potential is discretized as suggested in 



9] 



L, (j>n) - C = 



G(4> n ) - G(<t> n -!) 



where G' {(f) = ^j2V{<j>) - C, 



(21) 



then the last term of the Hamiltonian Eq. (|20|l reduces 
to (2/h)[G((f> n ) — G((f> n -i)] and it disappears in the tele- 
scopic summation. With the choice of Eq. (|2ip. the 
equations of motion derived from Eq. I|20() are 

d 

4>n = ~v(h, 0„_i, 0„)— — "0(/l, 4>n-l,4>n) 
0(f>n 

d 

-v(h, <j> n ,<f> n+ i)-Q—v(h, (f> n ,(f> n+1 ), (22) 

where, according to Eq. i|15|) and Eq. Q21[l. 

i-i G{cj) n ) - G(<j) n -x) 



i>(h, 4> n -i,<j) n ) 



•■"a Vn- 



(23) 



This lattice conserves the Hamiltonian (total energy) 

Obviously, static solutions to the lattice Eq. (1221 can be 
found from the two-point DFI v(h, <^ n _i, </>„) = and it 
is also clear that the Hamiltonian Eq. I|24l) is PNb-free. 

The energy-conserving lattice involves the integration 
constant C (through the function G) and this is different 
from what we had for the momentum-conserving model. 
The special case of C = yields the original model by 
Speight 9] , which supports a kink solution. Models with 
C =/= describe solutions different from the kink solution. 



Suppose that we use the same function V((/) n -i, (f> n ) 
in the DFIs Eq. {H} and Eq. (0 to construct 
different Di(h, <j) n -i, 4>n, 4>n+i) terms. Then equations 
U(h,4>n-i,4>n) = and v(h, <fi n -i, <j> n ) — are equiva- 
lent. Linear combination of those terms can be used to 
write the following PNb-free Klein-Gordon model 



^2biDi(h, 0„_i 



(25) 



where the constant coefficients satisfy the continuity con- 
straint 



5> 



1 when h — > 0, 



(26) 



assuming that bi can depend on h. Static solutions of 
Eq. (|25(l can be found iteratively from Eq. (|14f> or Eq. 

The model of Eq. H25fl can be generalized in a 
number of ways. For example, one can append terms 
which disappear in the continuum limit and ones that 
vanish upon substituting U(h, n -i> <Pn) — 0. Let us 
call such terms as O-terms. Furthermore, any term 
) can be modified by multiplying by 
a continuous function e(h, 4> n -i, 4>n, 4>n+i), which never 
vanishes and whose continuum limit is unity. Such mul- 
tiplication will not change either the continuum limit, or 
the static solutions of the model. 

Particularly we will study the PNb-free model of the 
form 



<t>n = e(h, <j> n ) 



U(h,<j> ) — U(h, 4> n -lAn) 



+P 



'Jn+1 



4>n)v{<t> n -l, <t>n 



(27) 



which is the lattice Eq. I|17fl modified by a multiplier 
e(h, 4>n) and augmented with the O-term with arbitrary 
coefficient p. This O-term was constructed to fit to the 
I3 invariant offered in ^| for the discretization of the <ft 4 
field. 

Finally we note that expressions similar to Eq. l|ltj|) 
can be derived from derivatives dv(x)/dx and dw(x)/dx 
and they can produce new PNb-free models. 

III. PNB-FREE DISCRETE ^ MODELS 

To give examples of Klein- Gordon lattices, we will dis- 
cretize the well-known field with the potential 



v{4>) = \d-4> 2 )\ 



(28) 



where the parameter A = ±1. The corresponding equa- 
tion of motion is 



+ Ai 



(29) 
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The above, so-called, 4> 4 equation supports moving pe- 
riodic solutions that can be expressed in terms of the 
Jacobi elliptic functions. The latter can be reduced to 
localized hyperbolic function solutions (when the elliptic 
modulus m = 1). Bounded solutions of this sort were 
discussed in the context of structural phase transitions 
|l4j and were used as the starting point for derivation of 
exact solutions to a discrete </> 4 model [lfj. In section 
| VIIII wc complete the list of the solutions presented in 
|l4Ul6j ] by providing also the set of unbounded solutions 
of the present model. 

The simplest discretization of Eq. I|29[l is of the form 



1 



t) n -l - 2(f) n + 4> n +l) + A I 



i>l), 



(30) 



and it possesses a PN barrier. However, as men- 
tioned above, one can construct PNb-free discrete 
Klein-Gordon models by discretizing the nonlinear term 
V'(<j)) typically over three neighboring points, V'(</>) — > 
V (4> n -i, 4> n , <f> n +i), in contrast to the classical discretiza- 
tion V'(<j>) — > V {<t> n )- The three-point discretization may 
be physically meaningful in some settings 17], but is also 
interesting from the more fundamental point of view of 
developing PNb-free discretizations and obtaining ana- 
lytically (or semi-analytically) explicit waveforms of their 
solutions. 

We discretize the potential of Eq. H28fl as follows 
A 



V(4>n-1>M = 1 (1- 



K-l4>n) ■ (31) 
Then, the DFIs, Eqs. (|14|) and (|15|l . become respectively, 



U(h, <f> n -l,<f>n) 



h 2 

^ (1 - 0„_i0„) 2 + C = 0, (32) 



v(h, <p n -i,(j> n ) 



^ (1 - K-i^n) 2 - C = 0. (33) 



It has been demonstrated that at C = Eq. I|32|) and 
Eq. support kink solutions HE1. For C ^ they 
support solutions different from kink. 
Let us denote 



A = Xh 2 , C = Ch 2 



(34) 



Substituting Eq. and Eq. (J33J into Eq. fgft with 
e(h, cj) n ) — 1 we arrive at the following PNb-free 4> A lattice 

4>n = A 2 0„ + \(f>n - ^<t>n(<t>n-l + 4>n+l) 



+ P{4>n+l - 4>n)\l ^ (1 - 4>n-l4>n) 2 - C 



P(4>n - 4>n-l)\l ^(l - <Pn(/>n+l) 2 - G. (35) 



Equation l|35(l gives a one-parameter (C) family of PNb- 
free models. For the special case of C = and A = 1 
we get one of the lattices derived in <j> n = (1 + 

h 2 p/V2)A 2 $ n + cj) n + {p/y/2 - l/2)^(^ n _i + (j) n +i) - 
p\/2(j) n -i(j) n (() n+ i, and this lattice supports the kink so- 
lution. Lattices with C / support solutions different 
from the kink. 

For p — we get from Eq. (fS^f> the C-independent 
model 



1 

+ Mn - 77^ 



n-1 
\ 2 

2 



20n + 4>n+l) 



(36) 



This non-Hamiltonian PNb-free 4 model conserves the 
momentum Eq. H18fl an d it will be referred to as the 
momentum-conserving (MC) model. 

Substituting Eq. into Eq. (23> wi th e(h, <j) n ) = 

1/(1 — Acf> 2 l /2) and p = we obtain another C- 
independent PNb-free model discovered in |lfij |. 



1 



20n + 4>n+l) + 



Ai 



1 - A02 /2 



(37) 



This model will be called the energy-conserving (EC) 
model because it possesses the Hamiltonian 



n 



+ 



2^ L r " ' h 2 

n 

where the potential V(<j> n ) is given by 



V{<t> n 



V(<t>n) = ~ 



2- A 
A 



In 



1 - 



HI 



(38) 



(39) 



In Fig. n]we plot the potential V(<j> n ) for h = 1 (i.e., 
A = A) in both cases A = 1 (solid line) and A = — 1 
(dashed line). For A < the potential is smooth and it 
has one minimum at <j) n = and two maxima at <f> n = ± 1 . 
For A > the potential has two minima at cj> n — ±1 and 
a maximum at <p n — 0; note that in the limit <f> n — > 
±\/2/A, the potential V(4> n ) — > +oo. 

It is not possible to plot an analog of Fig. ^ for the 
MC model since the relevant background forces are of 
many-body type (i.e., involve nearest neighbors) and are 
non-potential. 

PNb-free models given by Eq. (f53f) , Eq. lj36*|) , and Eq. 
(|37|l have exactly same static solutions derivable from 
DFI Eq. 



IV. EXACT STATIC SOLUTIONS FOR 
DISCRETE MODELS 

A. Solutions from nonlinear map 

To find all static solutions to the PNb-free models of 
Eq. ESJ, Eq. and Eq. (J27JI we solve the DFI of Eq. 
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FIG. 1: The on-site potential of the EC model of Eq. 
V((p n ), defined by Eq. (gU, for h = I and A = A = 1 (solid 
line) and A = A = — 1 (dashed line). For A < the potential is 
smooth and it has one minimum at <j> n = and two maxima 
at 4> n = ±1. For A > the potential has two m inima at 
(f>„ = ±1 and a maximum at (p n — 0; as (f>„ —> ±^/2/A the 
potential V(<j>„) — > +oo. 



ijniji: 



_ (2 - A)^ n _i ± VP 
P = 2A(1-0 2 „_ 1 ) 2 + 2C (A^_ x -2) 



(40) 



where <p n and n -i can be interchanged due to the sym- 
metry of the equation. Starting from any admissible "ini- 
tial" value (f> , by iterating Eq. ({4*0)1 and its counterpart 
written as <p n — l = f{,4 > n)i one can construct recurrently 
the static solution to both the MC model of Eq. 
and the EC model of Eq. 13 7|) . or to a linear combina- 
tion thereof. Arbitrariness in the choice of <f>Q implies 
the absence of PNb in these models, which has been also 
demonstrated in [Tfjl ]. 

As can be seen from Eq. I|4(JI) , once the values of C and 
A are fixed, there are certain restrictions on the choice of 
the values of <po. In particular, inadmissible initial values 
are those for w hich the denominator becomes zero, i.e., 
<f>Q =/= ±y2/A for A > 0. An exceptional case is that of 
A = 2, C = when an arbitrary sequence of ±1 is a 
solution of Eq. (J22J). Inadmissible values of </> are also 
ones for which T> < 0. The condition T> = leads to a 
biquadratic algebraic equation determining the borders 
of admissible region; the roots of this equation are: 




(41) 

Let us introduce the following notations for these roots: 

F 2 = -F 4 = x f(rth. (42) 



FIG. 2: Admissible region for the "initial" value <j>o in the 
nonlinear map Eq. 1401 for different values of C at A = 1. 
There are three inadmissible regions marked with "no" and 
two inadmissible lines (j>o / ±y2/A (horizontal dashed lines). 
On these lines the on-site potential V(4> n ) diverges (see Fig. 
0. Fi (i = 1, . . . , 4) designate different branches of the bor- 
ders of the admissible regions [see Eq. (011 and Eq. gg)]. 
Roots F\ and F2 merge at C = and at C — 4 — 8/ A. Roots 
F2 and F4 merge at C = A/2. Vertical dashed lines sepa- 
rate regions with different Jacobi elliptic function solutions. 
Vertical dotted line, situated (for the chosen parameters) at 
C* w —1.00, divides the region of the sndn/cn solution into 
two portions corresponding to two roots of the first equation 
in Eq. P . 



The admissible regions for the values of <j)Q of the non- 
linear map Eq. I|4(J|I are shown in Fig. [21 for different val- 
ues of C at A = 1 . The corresponding result for A = — 1 is 
presented in Fig. UJJ These graphs present a road map for 
constructing the various possible solutions of the above 
models. 

The symmetry of Eq. (|Tl2*|) suggests that the topology 
of the admissible regions is such that once started from 
an admissible value of 4>q , one cannot leave the admissible 
region iterating Eq. (|40|l , so that the static solution will 
surely be constructed for the whole chain. This is so 
because Eq. I|40() serves for calculating both back and 
forth points of the map, and if one is admissible, the 
other one is also admissible. 

Equation (|40|) possesses two roots, which means that 
for an admissible initial condition one can construct two 
different solutions, e.g., a kink and an antikink. When 
iterating, to keep moving along the same solution, one 
must take <j> n different from 4>n-i [if the roots of Eq. I|4(J|) 
are different]. Indeed, setting in the three-point static 
problem, Eq. (f> n -i — <t>n+ii wc find 



(2 - A)<t>„ 



(43) 



Fx 



-F a 



Comparing this with Eq. Q4Up. it is readily seen that 
the equilibrium in the three-point equation in the case 



7 




FIG. 3: Same as in Fig. Hbut for A = -1. Roots Fi and 
Fi merge at C = 0, and roots F2 and F4 merge at C = A/2. 
The dn solution, Eq. gSJ, is denned for A/2 < C < 0. The 
cn solution, Eq. 1481 . is defined for C < A/2, and this region 
is divided into two parts, (5\ and P2, each corresponding to a 
particular root of the first equation in Eq. 1481 . The border 
between these two regions is shown by the dotted line situated 
(for the chosen parameters) at C* ~ —2.96. 



and xq is the arbitrary initial position. Finally the inte- 
gers q, r, s specify a particular form of the solution. 

In the limit of m — 1, Eq. (|44l) reduces to the hyper- 
bolic function form 



±Atanh 9 (Z)cosh~ r ^(Z), 



(45) 



and, in the limit of to = 0, to the trigonometric function 
form 



±Asin«(Z)cos r (Z). 



(46) 



Substituting the ansatz Eq. l(4^|) into Eq. and 
equalizing the coefficients in front of similar terms we 
find that it can be satisfied for a limited number of com- 
binations of integer powers q,r,s. Solving the equations 
for the coefficients we relate the parameters of the solu- 
tion Eq. 144(1. (3 and A, to A and to and also find the 
relation between C, that enters Eq. and to. 

For some of the sets (q,r,s), e.g., for (q,r,s) — 
(1, —1, 0) and for (q, r, s) = (1, —1, —1), we obtain imag- 
inary amplitude A in the whole range of parameters. 

Essentially different, real amplitude solutions de- 
scribed by Eq. (|44|l . have the following form and are 
characterized by the following parameters: 

The sn solution, (q, r, s) = (1, 0, 0), 



of <f) n — 1 = (fin+i can be achieved only if T> — and the 
two roots coincide. As mentioned above, the latter re- 
quirement is equivalent to the condition that (fi n is on 
the border of the admissible region. 

As it can be deduced from Eq. ijlll) . the topology of 
the admissible regions presented in Fig. [21and Fig. [3]does 
not significantly change among different but positive and 
among different but negative values of A respectively; of 
course, one has to exclude the particular case of A = 
2 and also the continuum limit (see Sec. I Villi) . Here 
we will not discuss in detail the case of extremely high 
discreteness (A > 2), even though the analysis of this 
case docs not present any additional difficulties. 

To summarize, inside the admissible region, (C,(fio) 
(shown in Figs. [2Jand[2Jl, the nonlinear map of Eq. H4U|) 
generates static solutions to the MC model of Eq. 1|36|) 
and to the EC model of Eq. (|57|l . 



B. Jacobi elliptic function solutions 



A / 2m 

cn(/%)dn(/%) = l--, A=J— S n(f3h), 



2 V m 



The cn solution, (q, r, s) = (0,1, 0) 
cn(Ph) 



0<C<± (47) 



dn 2 (/3/i) 



A 



2m sn((3h) 



A dn(ph) ' 



C = A 



{I- A 2 



2 - AA 2 ' 
The dn solution, (q, r, s) — (0, 0, 1), 



n A 
-00 < C < — . 



dn(/3/i) 



cn 2 {ph) 



2^2 



C = A 



A aiiflhy 
(1 - A 2 ) 2 A 
2 - AA 2 ' 2 



The 1/sn solution, (q,r,s) = (—1,0,0), 



(48) 



-2 sn(/3fe) 
"cn(/3/i) : 

< C < 0. (49) 



Static solutions to the discrete PNb-free models of Eq. 
(|30ll and Eq. (|37|) have been reported 0] in the form of 
the Jacobi elliptic functions, sn, cn, and dn J2jj. Below, 
we will derive such solutions. 

The general form of the solutions is 

(fi n = ±Asn q (Z, 77i)cn r (Z, m)dn i (Z, to), 
Z = f3h(n + x ), (44) 

where < to < 1 is the modulus of the Jacobi elliptic 
functions, A and /3 are the parameters of the solution, 



A 12 
cn(f3h)dn(f3h) = 1 - -, A = J jsn(/3h), 



C=±{l-mA i ), 



< C < |. (50) 



The 1/cn solution, (q, r, s) = (0, -1, 0) 

cn(/3fe) 
dn 2 {fih) 



i _A A __/2i] ••,/,) s«.f iln 



A 



dn((3h) ' 



2 V 1 - m 



^<C*<oo. (51) 
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The sndn/cn solution, (q, r, s) = (1, —1, 1), 



men (Ph) + 1 — m 



en 2 (f3h) 



C 



A 
2 



A 

~r 

I -A 4 



A 



2 snM//'ldnM//'> 



A 



< C < 0. (52) 



It is worth making the following remarks: 

The solutions shown in Eq. (|4"5)l and Eq. have 
real amplitudes for A < while the others for A > 0. 

The solutions should be interpreted in the following 
form. For a given A, one can find j3 by solving the first 
equation in Eqs. 1)47152(1 . and then A from the second 
one. Substituting these values in Eq. (|44|l results in the 
static solutions of the original discrete model. 

The expressions for C in Eqs. 1(47152(1 link the elliptic 
Jacobi function solutions and the solution in the form 
of the nonlinear map, Eq. 1(40(1 . As for the other free 
parameter of the solutions Eq. H47I52|) . the arbitrary 
shift xq, its counterpart in the nonlinear map, Eq. ((401) . 
is effectively the initial value </>o- 

The solutions of Eqs. (|47I52II can be expressed in a 
number of other forms using the well-known identities 
for the Jacobi elliptic functions [23] • For example, shift- 
ing the argument by a quarter period, one can transform 
the sn solution to the form of cn/dn, or, applying the as- 
cending Landen transformation, to the form of sncn/dn. 
Mathematically, these three expressions look as different 
members of Eq. (|44|) . but physically they are indistin- 
guishable. 



V. ANALYSIS OF STATIC SOLUTIONS 

In what follows, we analyze the static solutions deriv- 
able from Eq. (|4"U|l discussing their relation to Eqs. l|T7|) - 
1(52(1. Since the topology of the admissible regions is dif- 
ferent for positive and negative A, these two cases are 
studied separately. 



A. A > case. 

In the sn solution, Eq. 1(47(1 , when the elliptic modulus 
increases from its smallest value to — to its largest 
value to — 1, the amplitude A increases from to 1, 
since A — F 2 [see Eq. (gTJ and Eq. ■ As a result, 
the parameter C monotonically decreases from A/2 to 
0. Thus, the sn solution is defined in the portion of the 
(C, </> )-plane, < C < A/2 and \<j) Q \ < F 2 (see Fig. EJ). 

The 1/sn solution, Eq. lf5U|) . is defined for < to < 1, 
and it is complementary to the sn solution since it is also 
valid for < C < A/2, but for \(f> \ > Fx (see Fig. 0. 

The 1/cn solution, Eq. (|51|l . is defined for < to < 1 
and it occupies the region A/2 < C < oo, |0 O | > F\ (see 
Fig. 0). 

The sndn/cn solution, Eq. 1(52(1 . is defined for unlim- 
ited (f> in the range 4-8/A<C<0 (see Fig. |2J). This 



solution is only valid for to* < to < 1, where, for fixed 
A, m*(h) is an increasing function of h and m*(0) = 1/2. 
For to < to* the first expression in Eq. 15211 does not 
have solutions for [3. For to* < to < 1, the equation 
has two roots, (3\ < @2- For the limiting value, to*, one 
can find the corresponding amplitude A* from the sec- 
ond expression of Eq. 1(52(1. and then C* , from the last 
expression. For the case of A = 1 presented in Fig. [21 we 
find to* w 0.933 and C* w -1.00. 

It should be noticed that we have not found a solution 
of the form of Eq. (glj valid in the range of C < 4 - 8/ A 
(portion marked with the question mark in Fig. It 
is likely that static solutions in this range cannot be ex- 
pressed in terms of the Jacobi elliptic functions because 
they do not survive in the continuum limit (see Sec. 
I Villi . However, the solution can easily be constructed 



from the nonlinear map in Eq. I|4U|I . 

Let us now discuss further several particular examples 
of the above solutions. First of all, in the limit m — ► 1 
[see Eq. (1451) ] . the sn solution, Eq. 1(4711 . reduces to the 
kink solution |16| . 



4> n = ± tanh[/3/i(n + Xq)], 



(53) 



while the 1/sn solution, Eq. 1)50(1 . to the solution called 
hereafter the "inverted" kink, 



±1 



tanh[/3/i(n + xq)] 



(54) 




In Eqs. ((53(1 and ((54(1 . x is the (arbitrary) position of 
the solution and tanh(/3/i) = \/ A/2. 

This limiting case corresponds to C = 0, for which 
a hcteroclinic connection is possible between the fixed 
points cj) n = — 1 and 0„ = 1 (see Fig. giving rise 
to the kink or the inverted kink. In this case, Eq. 1(40(1 
assumes the following simple form 



(55) 



where one can choose either the upper or the lower signs. 
The kink, Eq. (jSHJ, and the inverted kink, Eq. itBl^l . 
can be derived from Eq. ((55(1 taking initial values from 
1 0o | < 1 and 1 0o | > 1, respectively. 

In Fig. 0] we show (a) the kink and (b) the inverted 
kink solutions taking 4>o = \[2 — 1 and 0o = + 1 , 
respectively, for A = 1. 

In Fig. and for A = 1, we present two examples of 
solutions for positive and negative C close to 0. In par- 
ticular, for C = +2 x 10~ 5 , taking initial value 0o = 0, 
we obtain from the map of Eq. ((40(1 the solution pre- 
sented in panel (a). In fact, it is the sn solution close to 
the hyperbolic function limit having the form of a peri- 
odic chain of kinks and anti-kinks. On the other hand, 
for C = —6.75 x 10" 4 and <fio = 0, we obtain from the 
map of Eq. 1(40(1 the solution shown in panel (b). This 
is the sndn/cn solution close to the hyperbolic function 
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(a) 



(a) 
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FIG. 4: Solutions at C = 0: (a) kink and (b) inverted kink 
at A = 1. The solutions can be found from the nonlinear 
map Eq. 1551 with the initial conditions (j>o = v2 — 1 and 
4>o = v2 + 1 for (a) and (b), respectively. The kink and the 
inverted kink are also given by Eqs. lo.'-it and 1540 respectively, 
for xo — 0.5. Horizontal dashed lines, <f> n = ±y2/A, show 
positions of singular points of the potential V((j> n ) (see Fig. 

m 



FIG. 5: Solutions for C close to 0, at A = 1, obtained from the 
map of Eq. (a) <po = 0, C = 2 x 10" 5 and (b) 4>o = 0, 

C = —6.75 x 10 -4 . The solution in (a) is the sn solution, 
Eq. 1471 . and the solution in (b) is the sndn/cn solution, Eq. 
ll")2t . Horizontal dashed lines, </> n = ±\/2/A, show positions 
of singular points of the potential V{4> n ) (see Fig. 0. 



limit and has a form of a chain of kinks and inverted 
anti-kinks. 

In Fig. El and for A = 1, we present the 1/sn solu- 
tion, Eq. H50(l . (a) close to the hyperbolic and (b) at the 
trigonometric limits. The solution shown in panel (a) is 
a chain of inverted kinks and inverted anti-kinks. These 
solutions are obtained from the nonlinear map Eq. i|4U|) 
setting C = 1.23 xlO -9 , O = l + V2for (a), and C = 0.5, 
</-o = 2.45 for (b). 

In summary, we have shown that the well-established 
(hyperbolic and elliptic function) solutions of the model 
correspond to a very narrow region of the two-parameter 
(4> a ,C) admissible space. A natural question is what is 
the typical profile outcome stemming from the use of 
other pairs of <j>$ and C in Eq. (|4Ufl. Generically, upon 
testing the different regions of the admissible regime we 
have observed that arbitrary choices may lead to seem- 
ingly erratic solutions with very large amplitudes. A dif- 
ferent sign choice in the right hand side of Eq. (|4(JI) may, 
however, lead to a periodically locked tail structure. A 
simple example of such a solution is given below. 



B. 



A < case. 



and 



For negative A, we have the cn solution, Eq. 1)48 1). 
the dn solution, Eq. (|4*9l . 

Let us first start from the latter: When m increases 
from to 1 in Eq. I|49|) , the parameter C decreases mono- 
tonically from to A/2, i.e., this solution is defined in the 



'n n 
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FIG. 6: The 1/sn solution, Eq. (15011 . (a) close to the hy- 
perbolic and (b) at the trigonometric limits. We set A = 1. 
The solution in (a) consists of a chain of inverted kinks and 
inverted anti-kinks. To obtain these solutions from the non- 
linear map Eq. flU we set C = 1.23 x 10~ 9 , (f>o = 1 + V2 for 
(a), and C = 0.5, O = 2.45 for (b). 



"dn" portion of Fig. 01 

On the other hand, the cn solution of Eq. (|48|l is only 
valid for m* < m < 1, where m*(h) is an increasing 
function of h and m*(0) = 1/2. For m < m* the second 
expression of Eq. I|48|l does not have solutions for (3. For 
to* < to < 1, the equation has two roots, f3i < /?2- For 
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the limiting value, m*, one can find the corresponding 
amplitude A* from the second expression of Eq. H48f) and 
then C* from the last expression. When to increases from 
to* to 1 in Eq. Q48p. the parameter C of the nonlinear 
map Eq. (|40|l corresponding to the root f3i {fa) increases 
from C* (decreases from C*) to A/2 (to — oo). Thus, the 
cn solution occupies the rest of the admissible region for 
the case of A < (see Fig. [3Jl. In the case of A = — 1 
presented in Fig. we find m* » 0.873 and C* » -2.96. 

Both cn and dn solutions, in the limit m — > 1 (C — > 
A/2), reduce to a homoclinic to pulse solution (see also 
Figs, ^and |3 where this possibility is illustrated). This 
solution has the form: 

<J) n = ± A sech[/3/i(n + £ )], (56) 

where 



cosh(/3/i) = 1 - -, and A=d2--. (57) 

This is illustrated by Fig. where, taking A = — 1, we 
show the solution obtained from the map Eq. (|4UII at 
(a) C = A/2 - 10~ 8 (dn solution) and (b) C = A/2 + 
10~ 8 (cn solution), for initial value of 4>o = 10~ 4 . The 
figure clearly illustrates the two limits (to the left and 
to the right of C = A/2) and their correspondence to 
pairs of pulse-pulse solutions and ones of pulse-anti-pulse 
solutions, as one enters the two different regimes "dn" 
and "cn" of Fig. EJ 

C. Solutions with 4> n -i = 4> n +i 

Solutions of this special form can be expressed as cf> n = 
Acos(7rn) + B with constant A and B. For the sake of 
simplicity, we set B = and substitute the ansatz into 
Eq. (|36l) or Eq. (|37[1 to find the zigzag solution 



4> n = Acos(Trn), A=J—-1. (58) 

To obtain this solution from the map Eq. (|40|l . one has 
to set C = (A/2)(l + A 2 ) 2 -4A 2 //i 2 , which is the general 
condition for getting <fi n = — <^> n _i. Substituting here A 
from Eq. we get C = 4 (A - 2)/(Ah 2 ). For the case 

of A = 1 presented in Fig. [21 we have C = — 4. 

The zigzag solution is an exceptional one, as <f> n — 4>n-2 
for any n. However, as shown above, this is only possible 
when Eq. (|40l) has multiple roots, which means that the 
two-point static problem, Eq. is factorized. This 

solution does not have a counterpart in the continuum 
limit. 



VI. LINEAR STABILITY 

The static solutions of the discrete models Eq. (|36|) and 
Eq. (|37[1 are exactly the same, but the dynamical proper- 
ties of the two models are different. As the corresponding 
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FIG. 7: Solutions for A = -1 at (a) C = A/2 - 1(T 8 and (b) 
C = A/2 + 10~ 8 , for initial value of O = 10~ 4 . Since C is 
close to A/2 = —1/2, corresponding to m = 1, the dn solution 
in (a) and the cn solution in (b) look like chains of separated 
pulses given by Eq. 

stability analysis takes into account the dynamical form 
of each model, it will be carried out separately for each 
of the two. 

Let us first consider the MC model of Eq. (|36|l . In- 
troducing the ansatz 4> n (t) = </>° + s n (t) (where 0" is an 
equilibrium solution and s n (t) is a small perturbation), 
we linearize Eq. 1)36(1 with respect to e n and obtain the 
following equation: 

£n = T^{ e n—1 ~ + £n+l) + A£ n 

-^«) 2 (£n-l + £n+l) ~ A0M_ x + €+l)£n. (59) 

For the small-amplitude phonons, e n = exp(iKn + iu)t), 
with frequency oj and wave number k, around the uniform 
steady states (f>o = ±1 (A > 0), Eq. I|59|l is reduced to 
the following dispersion relation: 

- 2 = 2A+(^-2A)sin 2 (£), (60) 

while the spectrum of the vacuum solution for A < 0, 
< = 0, is 

^ Sin2 (i)- A - ( 61 ) 

For an arbitrary stationary solution </>° , stability is in- 
ferred in the MC model if the eigenvalue problem ob- 
tained from Eq. I|59l) by replacing e n with —uj 2 e n has only 
non-negative solutions ui 2 . Recalling that the MC model 
is a non-Hamiltonian one, the resulting eigenvalue prob- 
lem is a non-self- adjoint one involving a non-symmetric 
matrix. 
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Similarly, we obtain analogous expressions for the EC 
model of Eq. (|37[) . The linearized equation reads 



£n — J^( £ rL-i — 2e n + £ n +i) 

+ 2 ' A 2 + (A-6) W )' + A W 
[2-A(«)*f 

and the corresponding dispersion relation for the linear 
phonon modes around </>° = ±1 has the form: 



2 4A 4 2 
W = 2^A + 7? Sm 



(f 



(63) 



while that of vacuum solution </>„ = (for A < 0) is 



4 

7? 



• 2 / K 

sin - 

V2 



A. 



(64) 



The stationary solution (/>° is stable in the EC model of 
Eq. I|37fl if the self-adjoint eigenvalue problem obtained 
from Eq. (|62[1 by replacing e n with — Lu 2 s n has only non- 
negative solutions uj 2 . 

We note that all stable and unstable static solutions 
of the MC and the EC models, except for the zigzag so- 
lution Eq. (J5SJ, possess a zero-frequency mode. This is 
a consequence of the effective translational invariance of 
the discrete PNb-free models; this is related also to the 
freedom of selecting the free parameter xq in the corre- 
sponding solution expressions. 

Our aim here is not to study the whole bunch of so- 
lutions in the whole range of parameters, but rather to 
demonstrate the existence of stable solutions and also 
provide some examples where solutions, being stable in 
one model, may be unstable in another. The results 
of stability analysis for some characteristic solutions are 
summarized in Table [I] 

First, we consider the kink and inverted kink solutions 
given in Eq. (|53|l and Eq. i|54[l respectively: In a numeri- 
cal experiment, these solutions were placed in the middle 
of a chain of N = 200 particles with clamped boundary 
conditions. We have found that for the chosen set of pa- 
rameters, the kink is stable in both MC and EC models, 
while the inverted kink is stable in the EC and unstable 
in the MC model. 

In Fig. 121 we show the spectra of the kink and the 
inverted kink for their different positions with respect to 
the lattice, x$. In the figure, the horizontal lines show 
the borders of the phonon band of the vacuum cb n = ±1 
and the dots show the frequencies of the kink's internal 
modes lying outside of the band. In all cases, the spectra 
are shown as functions of xq, even though, admittedly, 
the kink in the MC model presented in (a) demonstrates 
very weak sensitivity of its spectrum to variations of Xq . 

Let us now consider the periodic solutions depicted 
in Fig. E{b) (sndn/cn) and in Fig. HJa) (1/sn). For 
these classes of solutions we used periodic boundary con- 
ditions and the length of the lattice was commensurate 
with the period of the solutions containing a number of 
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FIG. 8: Spectrum of (a) kink in MC model, (b) kink in EC 
model, and (c) inverted kink in EC model at C = 0, A = 1 
(h — 1). The inverted kink in MC model is unstable for the 
chosen parameters. Horizontal lines show the borders of the 
phonon bands. 



periods. Close to the hyperbolic function limit, these so- 
lutions are stable in the EC model but unstable in the 
MC model, following the stability of the building blocks 
(inverted kinks) associated with these structures. On the 
other hand, the 1/sn solution becomes unstable in both 
models at the trigonometric limit depicted in Fig. Elb). 
Nonetheless, we find it remarkable that some of these 
"apparently non-smooth" solutions (or maybe more cor- 
rectly, solutions apparently containing a non-smooth con- 
tinuum limit) may be potentially stable in the discrete 
setting. It would be wortwhile to analyze the stability of 
such structures in more detail and as a function of their 
elliptic modulus m. 

Next, we consider the zigzag solution, Eq. H58|) . which 
exists for < A < 4. It can be shown that this solu- 
tion can be stable in the EC model: Indeed, substituting 
Eq. (|58|l into the relevant eigenvalue problem, we ob- 
tain the dispersion relation for the linear phonon modes, 
lu 2 = (4//i 2 )sin 2 ( K /2) + 8/[h 2 (2 - A)]. The solution is 
stable when uj 2 are non-negative, i.e., when A < 2. Com- 
bining this with the existence condition we find that the 
zigzag solution exists and it is stable in the EC model 
for < A < 2. Note that this exceptional solution does 
not possess the zero frequency mode, as it can be seen 
from the above dispersion relation. The amplitude of the 
solution, ^/4/A — 1, is always greater than 1. The zigzag 
solution is always unstable in the MC model, a result 
that can similarly be demonstrated. 

The sn, cn, and dn solutions close to the hyperbolic 
limit were found to be generically unstable in both the 
MC and EC models. The instability of the pulse so- 
lution of Eq. (|56|l can be analytically justified in the 
present setting. It can be easily inferred from the invari- 
ance of the solution with respect to xo that the eigen- 
vector leading to the zero eigenvalue of Fig. |S] is propor- 
tional to d<j)/dx (where is the pulse solution profile). 
This is an anti-symmetric eigenvector, given the pulse's 
symmetric (around its center) nature. But then, from 
Sturm-Liouville theory for discrete operators [Tsj . there 
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should be an eigenvalue such that ui 2 < with a sym- 
metric (i.e., even) eigenvector, resulting in the instability 
of the relevant pulse solution. Our numerical simulations 
fully support this conclusion. The pulse solution was 
constructed by iterating Eq. gDJ f° r A < 0, C" = A/2, 
and < 4> < y/2 - A/2, or directly from Eq. JSBJ. We 
found that for A = — 1 the pulse solution is unstable in 
both models and over a wide range of the discreteness pa- 
rameter h, including the case of rather small h w 0.1. We 
have also confirmed that the instability mode is similar 
to the pulse profile (i.e., of even parity). 

Both the cn and dn solutions (see, e.g., Fig. 0) are 
also found to be unstable for the different parameter val- 
ues that we have used in both the MC and EC models. 
This is rather natural to expect given that their "building 
blocks" , namely the pulse-like solutions are each unstable 
(hence their concatenation will only add to the unstable 
eigendirections) . 

Finally, as far as the sn solutions are concerned, two 
basic instability modes were revealed. One resembles the 
instability of the pulse solution, when the whole struc- 
ture, being waved around the top of the background po- 
tential, tends to slide down to one of the potential wells. 
The other instability modes can be understood if one re- 
gards the sn solution as a set of alternating kinks and 
anti-kinks [see Fig. |Sfa)] interacting with each other by 
means of overlapping tails. Such a system is unstable 
because the kinks tend to annihilate with the anti-kinks. 

It is interesting to note that these stability results sug- 
gest a disparity with the earlier numerical findings of [lfj 
where the pulse-like solution was generically reported to 
be stable and analogous statements were made for the 
cn- and dn-type solutions. 

At this point, it is worth discussing a connection with 
NLS-type models mentioned in the introduction. In the 
case of such models, the aforementioned unstable eigendi- 
rcction (with an eigenvector similar in profile to the pulse 
itself) is "prohibited" by the additional conservation law 
of the I 2 norm, resulting in a separate neutral eigendi- 
rection with zero eigenfrequency. As a consequence, it is 
an interesting twist that the instability reported for such 
pulse-like solutions in scalar </> 4 models would be absent 
in their discrete NLS counterparts. 



VII. SLOW KINK DYNAMICS 

In |l5j . we have compared the spectra and the long- 
term dynamics of the kinks in the two PNb-free models, 
Eq. (|36|l and Eq. (|37|l . with that of the classical discrete 
model, Eq. I|30|l (see Fig. 1 of that paper). The case 
of very small kink velocities was not studied there. Here 
we would like to focus on this case to demonstrate the 
qualitatively different behavior for slowly moving kinks 
in the PNb-free and classical models. 

To boost the kink in the PNb-free models Eq. 
and Eq. (|37|l we used the following dynamical solu- 
tion corresponding to the multiple eigenvalue to 2 = 0, 



TABLE I: Results of stability analysis 



Solution 


MC model, 


EC model, 




Eq. TO 


Eq. EOT 


kink, Eq. J53J 


stable 


stable 


inverted kink, 


unstable 


stable 


Eq. El 






pulse, Eq. ll^t 


unstable 


unstable 


sn cn, dn 


unstable 


unstable 


close to hyper- 
bolic limit 






sndn/cn, 1/sn 


unstable 


stable 


close to hyper- 
bolic limit 







4> n (t) = + cte n , where </>° is a static kink solution, e n 
is the normalized translational eigenmode and c is the 
amplitude playing the role of kink velocity. This con- 
struction yields a more accurate approximate solution for 
very small c when linearized equations Eq. 159|) and Eq. 
(|62J) are accurate. Increasing c leads to a decrease in the 
accuracy of the dynamical solution used for boosting. On 
the other hand, we note in passing that there are tech- 
niques developed recently for Klein-Gordon 0, 0HJ and 
even nonlinear Schrodinger lattices that also allow 
the construction of finite speed, numerically exact trav- 
elling solutions in these classes of models. 

In Fig. [5] we present the lowest frequency normalized 
eigenmodes for the inter-site kinks in the three models, 
for the case of h = 1, A = 1 (A = 1). Both PNb-free mod- 
els, Eq. I|3tjfl an d Eq. (|37|) . have the same translational 
eigenmodes with oj — 0, shown by dots connected with 
solid lines. One can show by straightforward algebra that 
the PNb-free model modified by a non-vanishing multi- 
plier e(h,4> n -i,<f) n ,<f) n +i) as described in Sec. Ill Dl has 
the same translational eigenmode as the original model. 
Assuming that the two models are respectively: 

§n = f{4>n-l,4>n,4>n+l) (65) 
<j>n = e(4>n-l, <t>n, <t>n+l) X f(<t>n-l> <j>n, <l>n+l), (66) 

where the function e[x, y, z) ^ 0, then the steady state 
solutions 4>^i w iH satisfy f{(j^ l _ l , <j>^ , 4>n+i) — 0, for all 
n, while the corresponding linearization equations (using 
<t> n = 4>n + e n) ar e respectively of the form 



Of 



»n-l 

Of 



Of 



Of 



df_ 

0<t> n 



0(f>n+l 

Of 



e„+i, (67) 



0<Pn+l 



Hence, when solving the corresponding eigenvalue prob- 
lem (again substituting e„ by — w 2 e n ), for lo = 0, the 
eigenvalue problems become identical, both satisfying 



= 



Of 



0(j)n- 



Of 



Of 



'M+l 



£„+!, (69) 
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hence the coincidence of the corresponding eigenvectors. 
We are extremely grateful to an anonymous referee for 
this remark and its proof. We do note, in passing, also 
that while this justifies the coincidence of the zero eigen- 
frequency modes of the PNb-free models in Fig. it is 
also in tune with the results of Fig. |SI for nonzero eigen- 
frequencies to. The latter are not identical between the 
different models, as indicated by the left and middle pan- 
els of that figure. This is due to the differences between 
the corresponding eigenvalue problems of Eqs. H67[) an d 
, when uj^O. 

For the classical model, Eq. I|3U|I . the lowest-frequency 
mode has u> « 0.252, and it is shown in Fig. El by open 
circles and dashed lines. Actually, this mode is not a 
translational mode (since, strictly speaking, there is no 
translational invariance) but we use it to boost the kink. 
One can say that this will become the translational mode 
for this model in the continuum limit. 

We define the kink center of mass as 



S 



(70) 



The evolution of the kink coordinates is shown in Fig. 
1101 Kinks were boosted with two different amplitudes 
of the normalized lowest-frequency eigenmodes, c = 0.02 
and c = 0.08. Results for the PNb-free models practically 
coincide, as depicted by the dashed and solid lines for the 
MC and EC models, respectively. It is readily seen that 
kinks in the PNb-free lattices propagate with roughly 
constant velocities. 

The oscillatory trajectories in Fig. ^| correspond to 
the kink in the classical discrete 4> 4 model. The faster 
kink, boosted with c = 0.08, propagates along the lattice 
but its velocity gradually decreases. The kinetic energy 
of the translational motion is partly lost to the excitation 
of the kink internal mode with u ss 1.26, lying below the 
phonon frequency band. Higher harmonics of this mode 
interact with the phonon spectrum producing radiation, 
which also slows the kink down. An even more dramatic 
difference is observed for the slower kinks, boosted with 
c = 0.02. Here, the classical kink cannot overcome the 
PN barrier and can not propagate, oscillating near the 
stable inter-site configuration. On the other hand, the 
kinks in the PNb-free models are not trapped by the lat- 
tice and propagate due to the absence of the PN barrier. 
Alternatively, one can say that such waveforms can be 
accelerated by arbitrarily small external fields. 



.e 0.4 




FIG. 9: Lowest frequency normalized eigenmodes for inter- 
site kinks. Both PNb-free models, Eq. (O and Eq. 157) . 
have the same translational eigenmode with ui = 0, shown 
by dots connected with solid lines (see also the relevant dis- 
cussion in the text). For the classical model, Eq. 13011 . the 
lowest-frequency mode has uj m 0.252 (shown by open circles 
and dashed lines). Results for h — 1, A = 1 (A = 1). 




VIII. SOLUTIONS FOR CONTINUUM 4 FIELD 

In the continuum limit, h — ► 0, the borders of the 
admissible region, Eq. lBT|) . become 

(0g) 1(2 = l±^?. (71) 

In Fig. ^]we plot the admissible regions for (a) A = 1 
and (b) A = — 1 . The topology of the admissible regions 



FIG. 10: The kink coordinate, S, as a function of time for 
kinks propagating in the PNb-free MC (dashed line) and EC 
(solid line) models, as well as in the classical discrete <j> 4 model 
(oscillatory line). The kinks were boosted with two different 
amplitudes of the normalized lowest-frequency eigenmodes, 
c = 0.02 and c = 0.08. The faster classical kink is able 
to propagate, while the slower one is not, since it cannot 
overcome the PN barrier. Kinks in PNb-free models are not 
trapped by the lattice and can therefore propagate. 
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FIG. 11: Admissible regions for the 4> 4 field for (a) A = +1 
and (b) A = — 1 obtained as the continuum limits (h — > 0) of 
those presented in Fig. [5] and Fig. [3] respectively. In each 
panel there is one inadmissible region marked with "no". 



for the 4 field is simpler than the one pertaining to the 
discrete models. In the continuum limit there exists a 
sole inadmissible region for both cases A = ±1, while 
three inadmissible regions exist for the discrete models 
in the case A > 0. One more simplification is that the 
domains of the sndn/cn and cn solutions do not split 
into two parts since the smaller root /3i disappears in the 
continuum limit. Particularly we note that the region 
marked with the question mark in Fig. |3 C < 4 — 8/ A, 
disappears in the continuum limit. That might be the 
reason why we failed to find a Jacobi elliptic function 
expression for the discrete static solutions in this case. 



Static solutions Eqs. 1(47152(1 obtained for the discrete 
models have their continuum counterparts as the travel- 
ing solutions to the 4 field, Eq. 1)29(1 . The general form 
of the solutions is 



4>(x,t) = ±Asn q (z, m)cn r (z, m)dn s (z, to), 
z = P (x + Xq — ct) , 



(72) 



where < to < 1 is the modulus of the Jacobi elliptic 
functions, A and (3 are the parameters of the solution, 
xq is the arbitrary initial position and < c < 1 is the 
velocity of the solution. The integers q, r, s once again 
specify the particular form of the solution. 

Continuum analogues of Eqs. 1(47(1 - 1(52(1 have the fol- 
lowing form and are characterized by the following pa- 
rameters: 



The sn solution, (q, r, s) 

= 




A 



to)(1 



2 \ m 



(1,0,0), 



2m 



0<C<-. 



The cn solution, (q, r, s) 

(3 = 



-A 



(2m- l)(l - c 2 ) 



(0,1,0), 
, A = 



2m 



2m- 1 



C 



„ A 
-oo < C < —. 



The dn solution, (q, r, s) 

= 



-A 



(0,0,1), 
. A = 



(2-to)(1-c 2 )' " V 2- m 
A 2 ) 2 

The 1/sn solution, (q,r,s) — (—1,0,0) 



±<C<0. 



A 



(1 + to)(1-c 2 )' 



A = 



1 



C=±{l-mA% 0<C<^. 



The 1/cn solution, (q,r,s) = (0,-1,0) 



A 



C 



(1 - 2m) (1 - c 2 )' 

A A mA 4 
2 



A 



2(1 



1 - 2m 



1 



A 



1 — to J 2 
The sndn/cn solution, (g, r, s) = (1, —1, 1), 

1 








A 



V2to - 1 
-oo < C < 0. 



(73) 



(74) 



(75) 



(76) 



< C < oo. (77) 



(78) 



The above six solutions can be rewritten in a great 
variety of forms using the properties of the Jacobi ellip- 
tic functions [2(|. However, we believe that they may 
constitute the full list of physically different solutions to 
the continuum 4> equation since they fill the whole two- 
parameter space, (C, </>), obtained as the continuum limit 
of corresponding space of the discrete models. 

All solutions are conveniently parameterized by a sin- 
gle parameter — oo < C < oo for A > 0, and — oo < C < 
for A < 0, as it is presented in Fig. ITT1 

The solutions in Eqs. ((73|) - ((75(l are bounded while the 
other ones are unbounded. The solutions in Eq. I(74|) 
and Eq. ((75(1 are defined for A < while the others for 
A > 0. The solutions in Eq. (JHJ and Eq. O are valid 
for 1/2 < to < 1, the solution in Eq. 1(77(1 is valid for 
< to < 1/2, while the other solutions for < to < 1. 
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IX. CONCLUSIONS 

In the present paper we have shown that the reduc- 
tion of the static problem of a discrete Klein-Gordon 
(and by extension of the standing wave problem of a dis- 
crete NLS) equation to a two-point problem is a power- 
ful tool for obtaining all possible static solutions of the 
corresponding model. We have applied this general idea 
to a momentum conserving and an energy conserv- 
ing 0| discretization of the <f> 4 model analyzing the full 
two-parameter plane of solutions and giving a natural 
parametrization for it. In particular, we have examined 
the admissible regions of the field value at a given point 
and of the constant entering the two-point function per- 
tinent to the model. We have specifically illustrated how 
to use different choices of these two parameters to ob- 
tain not only the well-known hyperbolic function solu- 
tions and the established elliptic function solutions, but 
also novel classes of solutions including the inverted (non- 
monotonic) kink solutions presented herein and multi- 
kink generalizations thereof. We performed such com- 
putations both for the attractive (focusing) and for the 
repulsive (defocusing) types of nonlinearity. 

The presented methodology has the significant advan- 
tage over earlier work that by introducing the integration 
constant in the discretization of the first integral of the 
static part of the PDE, it allows one to construct the full 
family of the static solutions of the corresponding model. 
Earlier work had implicitly set this additional free pa- 
rameter to 0; this choice was sufficient for obtaining the 
important hyperbolic function solutions, but the present 
formalism systematically illustrates how to generalize the 
latter. 

We have also derived the continuum analogues of the 
discrete solutions and found that they fill the whole space 



of parameters obtained in the limit of h — > from the 
parameter space of the discrete models. 

In addition, we have examined some of the key stability 
features of the solutions obtained in the various models, 
illustrating the different stability properties of the models 
considered herein (even if their static solutions are iden- 
tical). We have obtained interesting stability properties, 
including, e.g., the counter-intuitive stability of the in- 
verted kink and of some of its periodic generalizations 
for the EC model. 

The present study generates a variety of interesting 
questions. For instance, it would be relevant to examine 
whether the solutions presented herein (including the in- 
verted ones) have counterparts in the "standard" discrete 
</> 4 model and to analyze their respective stability. It 
would also be relevant to examine how the stability (and 
existence -since some of them may not survive the contin- 
uum limit-) of such solutions depend on the lattice spac- 
ing h. Finally, these and related (e.g. stabilit y) ques tions 
in the context of discrete NLS lattices (see [23, Hj, |25j) 
would also be equally or even more (given the multitude 
of relevant applications of the latter model) interesting 
to answer. Such studies are currently in progress and will 
be reported elsewhere. 
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